Superconductivity in disordered thin films: giant mesoscopic fluctuations 
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We discuss intrinsic inhomogeneities of superconductive properties of uniformly disordered thin 
films with large dimensionless conductance g. It is shown that mesoscopic fluctuations, which 
usually contain a small factor 1/g, are crucially enhanced near the critical conductance g C F 2> 1 
where superconductivity is destroyed at T = due to Coulomb suppression of the Cooper attraction. 
This leads to strong spatial fluctuations of the local transition temperature and thus to percolative 
nature of the thermal superconductive transition. 
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Since the very early stage of superconductivity the- 
ory it is known [lj, 2] that the superconducting transition 
temperature, T c , is insensitive to the rate r _1 of elastic 
impurity scattering, i.e. it does not depend on the pa- 
rameter tT c /H. This statement known as the "Anderson 
theorem" is valid provided that both (i) Coulomb inter- 
action effects and (ii) mesoscopic fluctuations are negli- 
gible. However, in sufficiently disordered metals, close to 
the Anderson localization transition, these effects become 
important and the Anderson theorem is violated. 

In disordered samples, Coulomb repulsion enhanced 
due to diffusive character of electron motion |3| leads to 
the suppression of T c with the increase of disorder (cf. 
Ref. 0for a review). For two-dimensional (2D) thin films 
with the dimensionless conductance g = 2irh/e 2 R\j 3> 1, 
the first-order perturbative correction had been calcu- 
lated in Refs.0, and the general expression for T c (g) was 
obtained by Finkelstein 6]. In the leading order over 
g -1 / 2 his result reads: 
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where T c q is the transition temperature in the clean 
(g — ► oo) system, and r* = max{r, r(d/l) 2 }, with d being 
the film thickness and I = vpr being the mean free path. 
According to Eq. £[J, T c vanishes at the critical conduc- 
tance g C F = ln 2 (S/T c or*)/(27r) (which is supposed to be 
large enough for the theory to be self-consistent). 

The Finkelstein's theory is an extended version of 
the mean-field theory of superconductivity which takes 
into account that the effective Cooper attraction ac- 
quires (due to Coulomb interaction and slow diffusion) 
a negative energy-dependent contribution. Within this 
( "fermionic" ) mechanism, vanishing of T c is accompanied 
by vanishing of the amplitude of the superconductive or- 
der parameter A. Another ("bosonic") mechanism of T c 
suppression is due to phase fluctuations of the order 
parameter. This mechanism seems to be adequate mainly 
for structurally inhomogencous superconductors (granu- 
lar films of artificial arrays) with well-defined supercon- 
ductive grains interconnected by weak links. Below we 



will see, however, that phase fluctuations inevitably be- 
come relevant for homogeneously disordered films with 
the conductance g close to its critical value g c p. 

Another phenomenon important for 2D conductors is 
known as mesoscopic fluctuations 8] and is due to non- 
local interference of electron waves scattering on impuri- 
ties. It was recognized by Spivak and Zhou |9j that sim- 
ilar fluctuations are pertinent also for the Cooper paring 
susceptibility, K(r, r') = (K (r — r')) + 6K(r,r'), which 
enters the BCS self-consistent equation 



A(r) = ^J K(r,r')A(r')dr', 
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where A is the dimensionless Cooper coupling constant 
and v is the single-particle density of states per spin. 
Equation with the exact disorder-dependent kernel 
K(r, r') possesses localized solutions for A(r) above the 
mean-field transition line. They describe droplets of su- 
perconducting phase which nucleate prior to the transi- 
tion of the whole system. Since the relative magnitude 
of mesoscopic fluctuations of SK(r, r') is of the order of 
1/g and is small for a good metal, the effect of localized 
droplets on the zero-field superconductive transition is 
negligible and the transition width is determined by ther- 
mal fluctuations. Contrary, at low temperatures near the 
upper critical field H C 2(0), thermal degrees of freedom are 
frozen out and mesoscopic fluctuations are fully respon- 
sible for the width of the field-driven superconductor- 
normal-metal (SN) transition 0,^3- Still, the relative 
magnitude of H C 2 (0) shift and of the transition width due 
to mesoscopic fluctuations are of the order of 1/g <C 1 as 
long as Coulomb effects arc neglected. 

The goal of the present Letter is to develop a combined 
theory of the superconductive transition in 2D disordered 
films, which takes into account both Coulomb effects and 
mesoscopic fluctuations. Our main result is the expres- 
sion for the relative smearing 5d = ST C /T C of the zero-field 
transition due to formation of localized islands: 
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where ad « 0.4. In the vicinity of the Finkelstein's critical 
point, at g — g C F ^5 1, formation of islands dominates over 
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FIG. 1: Diagrams for the GL free energy: (a) |A| 2 term; (b) 
|A| 4 term; its central part is the Hikami box. 
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FIG. 2: Equation for the cooperon in the fluctuating electric 
field. The Coulomb vertex C is given by Eq. @. The kernel 
(K ek (r — r')) is obtained after summation over e m . 



thermal fluctuations characterized by the Ginzburg num- 
ber Gi= (n/8g) |H|. Moreover, in the very close vicinity 
of the quantum critical point, at g — g C F ^ ^/gcF, fluctu- 
ations of the "local transition temperature" become large 
on the absolute scale, ST C ~ T c , and the superconductive 
state becomes strongly inhomogeneous in the absence of 
any pre- determined structural granularity [l2j |. 

We emphasize that mesoscopic fluctuations are mini- 
mal fluctuations which are inevitably present in any dis- 
ordered system. In real samples, their effect may be en- 
hanced by various types of structural inhomogeneities. 

Ginzburg- Landau expansion. — We begin with deriving 
the Ginzburg-Landau (GL) expansion in the vicinity of 
the Finkclstein transition temperature 0). The GL free 
energy for the static order parameter has the form: 
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where the first term is the disorder-averaged contribu- 
tion, and the last term accounts for mesoscopic fluctua- 
tions; its form will be found later. The diagrams for the 
disorder-averaged free energy |l3| are shown in Fig. ^ 
Apart from the standard impurity averaging, cooperons 
should be averaged over fluctuations of the electric field 
(this is shown by the gray rectangle). The order param- 
eter always enters in the combination (K e (r — r'))A(r'), 
where we introduced the reduced Cooper kernel: 



(^(r-r'))=T^(G efc , em (r,r')G_ 
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with G £k y£m (r, r') being the Matsubara exact Green func- 
tion which, in the presence of the fluctuating electric field, 
depends on two energy arguments. After summation over 
energy Eq. (J3J) gives the (averaged) pairing susceptibility: 
(if(r-r'))=T£ £ (K E (r-r')>. 

The kernel (K e (r— r')) obeys the linear equation shown 
diagrammatically in Fig. El T° write it in a compact form 



we define the cooperon screening factor w q (e): 
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which shows how the free metallic cooperon gets mod- 
ified by the fluctuating electric field (D is the diffusion 
coefficient). The screening factor satisfies the equation 
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(7) 



D < ?+2\e n \' 
where (for thick films with d > I replace t by t,) 
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is the disorder-enhanced Coulomb interaction vertex [LH 
proportional to the return probability, with 6{x) being a 
unit step function. 

Equation (JJJ) looks pretty similar to the equation for 
the energy-dependent Cooper vertex r e)E ' considered in 
Ref.Q 

It can be easily solved with logarithmic accuracy 
[l4j |. Introducing a new variable ( = ln(l/er), one finds 
for the zero- momentum limit of the screening factor : 



w(() = w {e) — cosh(AgC) — tanh(A 9 (T) sinh(A s C), (9) 

where = ln(l/Tr), and X g = l/^/2ng is the Finkel- 
stein's fixed point. The function wo(T) decreases from 1 
at e ~ 1/r down to w((t) = 1/ cos1i(A 9 Ct) at e ~ T. 
The coefficients in the GL free energy J2J are given by: 
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where O = 7C(3)^/(87r 2 T c 2 ) and 7 = nvD/(8T c )axe the 
standard coefficients for dirty superconductors |l3j . and 
A* is the running Cooper coupling constant at the en- 
ergy scale r _1 . The Matsubara sums in Eqs. (pb) and 
(ll(Jc|) converge at the thermal scale. Therefore the coef- 
ficients (3 and 7 contain the screening factors evaluated 
at Qr- Contrary, the coefficient a is determined by all 
energies e < 1/r. The integral in Eq. (|10a|l is given by 
A" 1 tanh(A g (r), and solving a(T c ) = one immediately 
recovers the Finkel' stein expression (|TJ for T c . Taking 
the derivative near T c we find a = i/(T /T c — 1)w 2 (£t)- 

Thus, disregarding mesoscopic fluctuations, we see 
that the A field always enters the GL expansion in 
the combination with the screening factor w(Ct) = 
1/ cosh(A ff (r). If we define A = Aw((t c ) then the GL 
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FIG. 3: One-loop diagrams for the 4-cooperon/diffuson col- 
lision vertex M Ei (q, q'). The diagrams (a) and (c) have their 
symmetric counterparts. 



expansion for A will acquire the standard form with the 
coefficients a = v(T/T c — 1), 0o and 70. As a corollary, 
the Ginzburg number appears to be unaffected by the 
Coulomb repulsion: Gi = n/(8g). 

Quasiparticle spectrum. — Incorporating the Coulomb 
screening factors Wq(e) into the standard Green func- 
tion formalism, one finds that a quasiparticle propagat- 
ing with the energy e feels an effective pairing potential 
A e ff(£) = w(()A. The gap in the spectrum, E gap , should 
be found as a solution to E gap = iy[ln(l/^g ap r)]A. Deep 
in the superconducting phase, at T c — T ~ T c , it coin- 
cides with A. On the other hand, excited particles with 
£ > ^gap see a larger value of A e ff(e), and high-energy 
particles with e > 1/t feel the bare value of A. 

In the regime of strong Coulomb suppression of super- 
conductivity, when T c -C T c q and w(Ct c ) <C 1, the bare 
A significantly exceeds the screened A. This enhance- 
ment of A was irrelevant for the GL expansion since the 
sums over Matsubara energies which determined the GL 
coefficients converged at the thermal scale. We will see 
below that this is not the case for mesoscopic fluctuations 
where the bare value of A comes into play. 

Mesoscopic fluctuations of the pairing susceptibility. — 
In order to calculate the correlation function of the 
pairing susceptibility, one has to draw two diagrams 
for K(r, r') (see Fig. QJi) and connect their diffusive 
modes by impurity lines. In general, the variance 
(SK(ri,r2)5K(r 3 ,r4)} is a complicated function of — 
Tj decaying at the scale of the thermal length Lt — 
\J D I (2nT). On the other hand, close to T c the order 
parameter varies at the scale of the coherence length 
£(T) = L Ty /T c /{T -T c ) > Lt. Therefore, in the vicin- 
ity of the superconducting transition the fluctuations of 
K(r, r') are effectively short-ranged and characterized by 
the single number 



C = J (6K(r 1 ,r 2 )SK(r 3 ,r 4 )) dr 2 dr 3 dr 4 
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Taking into account all possible correlations between dif- 
fusive modes in K(r,r') we find 

C = T " E ( LI ^) ^ M ei (q, q'). (12) 
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where IZ^ is an operator acting on an arbitrary function 
X(q) as 
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Eij = £i + Ej, and M Ei (q, q') is the 4-cooperon/diffuson 
collision vertex shown in Fig. |3| (with the proper construc- 
tion of internal Hikami boxes by drawing additional im- 
purity lines being implied) . The internal diffusive modes 
can be diffusons or cooperons. In the vicinity of H C 2(0), 
the diagram (b) was considered in Refs. 0, and the dia- 
gram (c) was analyzed in Ref. ^| At zero magnetic field 
the vertex is calculated elsewhere 1 151: 



M £i (q,q') = 
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The first term in Eq. (|13fl refers to the cooperons in 
the ladder Q shown in Fig. |3 while its second term 
refers to the cooperons/diffusons which are responsible 
for the return probability in the vertex C given by Eq. ©. 
Some summations over energies in Eq. (|12(l saturate at 
the thermal scale, whereas summations associated with 
return probability are logarithmic extending up to the 
high-energy cutoff r _1 : 
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The second term in the square brackets in Eq. 1)15(1 is due 
to mesoscopic fluctuations of return probability. Calcu- 
lating the sum as an integral over £ with the help of 
Eq. I© we find that this term is equal to 1 — wq(T), 
so that the total expression in the square brackets in 
Eq. (|15|l is equal to 1. Thus, in the regime of strong 
Coulomb suppression of superconductivity, the pairing 
susceptibility fluctuates mainly due to mesoscopic fluc- 
tuations of return probability in the vertex (JSJ . Writing 
C = C w^T), we get 
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Superconductor with fluctuating T c . — Short-range 
mesoscopic fluctuations of (5if(ri, r 2 ) are equivalent to 
local fluctuations of the transition temperature ST c (r) 
which can be described by the following term in the GL 
free energy: 



= / <55(r)|A(r)| 2 dr, 
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where (Sa(r)Sa(r')) = CoS(r—r'). Superconductors with 
local fluctuations of T c were considered previously within 
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the phenomenological approach by Ioffe and Larkin [17j , 
where the three-dimensional case was mainly discussed. 
Generalizing their results to the 2D case, we find that 
the relative smearing of the superconductive transition 
due to frozen-in mesoscopic fluctuations is given by Sd — 
Co/[12T c 7o(ciao/dT)] (note that in the 2D case the nu- 
merical coefficient in the exponent of Eq. (29) in Ref. fl7l 
is equal to 11.8, cf. 18]). Taking C at T = T c and us- 
ing Eq. 0) we obtain a surprisingly simple expression (jjJJ 
with a d = 28C(3)/3tt 3 . 

The increase of Sd near the critical conductance g C F can 
be understood in terms of renormalization of the Cooper 
attraction constant A. At low energies, the latter acquires 
a negative Coulomb contribution proportional to the re- 
turn probablity ~ g~ x ln(l/er). Mesoscopic fluctuations 
of g lead then to fluctuations of A, whose relative effect 
grows with decreasing A: Sd ~ ST C /T C — SX/X 2 . 

Equation J2J predicts that for g—g C F ^ 1, the disorder- 
induced broadening of the transition dominates over the 
thermal width: Sd > Gi. In such a situation, the macro- 
scopic superconductive transition occurs via formation of 
small superconductive islands of size Lp = yD/ (2irT c ) 
surrounded by the normal metal state. With the tem- 
perature decrease, the density of these islands and the 
proximity-induced coupling between them grows until 
a percolation-type superconductive transition ^3] takes 
place. At sufficiently low temperatures, T < T c (l — Sd), 
the superconductive state becomes approximately uni- 
form, with weak spacial variations in the amplitude of 
the order parameter |A|. 

Another situation occurs in the closest vicinity of the 
critical conductance, g — g c p < 1/ g c p: here local fluctu- 
ations of T c are large on the absolute scale, and strong 
inhomogeneity of the superconductive order parameter 
persists down to T <C T c . As a result, both thermal and 
quantum fluctuations of phases of superconductive or- 
der parameters on different superconductive islands are 
strongly increased. In other terms, in the close vicinity 
of the critical conductance g C F, the bosonic mechanism 
of superconductivity suppression becomes relevant. 

Inhomogeneous distribution of |A(r)| is known [Tgt l20| 
to smear the gap in the excitation spectrum of a super- 
conductor. We expect this effect to be very strong for 
9 ~ 9cF- 

The result Q indicates that str ong en hancement (in 
comparison with the results of Refs.loLllOl) of mesoscopic 
fluctuation effects in H C 2 behavior at low temperature 
should be expected at g » g c p. This problem needs 
special treatment since correlation length of mesoscopic 
fluctuations diverges at T — > 0, and short-range approxi- 
mation employed for the determination of C in Eq. i|ll|) 
becomes inappropriate [l6j. It is quite clear however, 
that long-range features of mesoscopic disorder at T — > 
can only increase the effective width of the field-driven 
T = SN transition in comparison with the width Sd of 
the zero-field transition driven by temperature. An ex- 



tension of the present approach to magnetic-field-induced 
transition near H C 2(0) will be considered separately. 

To conclude, we demonstrated that strong inhomo- 
geneities of superconductive state can be induced by rel- 
atively weak (~ <? _1 ) mesoscopic fluctuations, which lead 
to spatial fluctuations of the effective Cooper attraction 
constant. As a result, a nominally uniformly disordered 
film may appear as a granular one in terms of its super- 
conductive properties. 
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